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Two major themes in the physics of condensed matter are quantum 
critical phenomena and unconventional superconductivity. These 
usually occur in the context of competing interactions in systems 
of strongly-correlated electrons. All this interesting physics comes 
together in the behavior of the recently discovered iron pnictide 
compounds that have generated enormous interest because of their 
moderately high-temperature superconductivity. The ubiquity of an 
tiferromagnetic ordering in their phase diagrams naturally raises the 
question of the relevance of magnetic quantum criticality, but the an- 
swer remains uncertain both theoretically and experimentally. Here 
we show that the undoped iron pnictides feature a novel type of 
magnetic quantum critical point, which results from a competition 
between electronic localization and itinerancy. Our theory provides 
a mechanism to understand the experimentally-observed variation of 
the ordered moment among the undoped iron pnictides. We sug- 
gest P substitution for As in the undoped iron pnictides as a means 
to access this new example of magnetic quantum criticality in an 
unmasked fashion. Our findings point to the iron pnictides as a 
much-needed new setting for quantum criticality, one that offers a 
new set of control parameters. 
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approach is motivated by the phenomenological and the- 
oretical evidence that the parent iron pnictide is a "bad 
metal" [9, 17, 18]. Accordingly, we formulate our consider- 
ations in terms of an incipient Mott insulator: the electron- 
electron interactions lie close to, but do not exceed the critical 
value for the insulating state. Within this picture, the elec- 
tronic excitations comprise an incoherent part away from the 
Fermi energy, and a coherent part in its vicinity. The incoher- 
ent electronic excitations are described in terms of localized 
Fe magnetic moments, with frustrating superexchange inter- 
actions. The latter have been discussed earlier by two of us [9] 
and others [16]. This division of the electron spectrum is a 
simple and convenient way of analyzing the complex behavior 
of a bad metal close to the Mott transition, whose spectrum 
exhibits incipient upper and lower Hubbard bands and a co- 
herent quasiparticle peak at the Fermi energy [19]. 

The coupling of the local moments to the coherent elec- 
tronic excitations competes against the magnetic ordering. 
A magnetic quantum critical point arises when the spectral 
weight of the coherent electronic excitations is increased to 
some threshold value. 



Quantum criticality in the pnictides 
The recent discovery of copper-free high-T c supercon- 
ductors has triggered intense interest in the homologous iron 
pnictides. The parent compound of the lanthanum-iron ox- 
yarsenide, LaOFeAs [1], exhibits a tetragonal-orthorhombic 
structural transition and long-range antiferromagnetic or- 
der [2] . Electron doping, via fluorine substitution for oxygen, 
suppresses both and induces superconductivity. Other fami- 
lies of the arsenide compounds show a similar interplay among 
structure, antiferromagnetism and superconductivity. These 
include the oxyarsenide systems obtained through replacing 
lanthanum by other rare-earth elements such as Ce, Pr, Nd, 
Sm, and Gd [3, 4, 5, 6], as well as oxygen-free arsenides, such 
as BaFe 2 As 2 [7] and SrFe 2 As 2 [8]. 

The existence of the antiferromagnetic state naturally 
raises the possibility of carrier-doping-induced quantum phase 
transitions in the iron pnictides [9, 10, 11], but the situation 
is not yet certain. Theoretically, the evolution of the Fermi 
surface as a function of carrier doping is not yet well under- 
stood, and this limits the study of quantum criticality. Ex- 
perimentally, earlier measurements in LaOi_ x F x FeAs [1] and 
SmOi_ x F x FeAs [12] show a moderate suppression of the mag- 
netic/structural transition temperature(s) as x is increased; 
beyond x of about ~ 7%, the transitions are interrupted by 
superconductivity. Further experiments have led to conflicting 
reports for the first-order or second-order nature of the carrier- 
induced zero-temperature magnetic and structural phase tran- 
sitions [13, 14, 15]. 

We propose that an alternative to a possible doping- 
induced quantum phase transition is one that is accessed 
by changing the relative strength of electron-electron corre- 
lations. Thus we suggest that the iron pnictides may ex- 
hibit a new example and setting for quantum criticality. Our 



The electron spectrum 

The incoherent and coherent parts of the single-electron spec- 
tral function are illustrated schematically in Fig. 1. The cen- 
tral peak describes the coherent itinerant carriers; these are 
the electronic excitations that are responsible for a Drude op- 
tical response and that are adiabatically connected to their 
non-interacting counterparts. The side peaks describe the in- 
coherent excitations, vestiges of the lower and upper Hubbard 
bands associated with a Mott insulator that would arise if 
the electron-electron interactions were larger than the Mott- 
localization threshold. Each of the three peaks may in general 
have a complex structure due to the multi-orbital nature of 
the iron pnictides. The decomposition of the electronic spec- 
tral weight into coherent and incoherent parts is natural for a 
metal near a Mott transition [20, 19]. 

We use w to denote the percentage of the spectral weight 
lying in the coherent part of the spectrum. A relatively small 
w may be inferred for the iron pnictides, because the Drude 
weight seen in the optical conductivity [21, 22, 23] is very 
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small (on the order of 5% of the total spectral weight inte- 
grated to about 2 eV). A small w corresponds to an interac- 
tion strength sufficiently large that the system is close to the 
Mott transition, albeit on the metallic side; this implies a large 
electron-electron scattering rate, consistent with the observed 
large electrical resistivity (on the order of 0.5m£7 ■ cm for single 
crystals and 5mf2 • cm for polycrystals) at room temperature. 
In terms of electrical conduction, the iron pnictides are sim- 
ilar to, e.g., V2O3, a bad metal (with a room temperature 
resistivity [24] of about 0.5 mS7 ■ cm) that is known to be on 
the verge of a Mott transition, and is very different from, e.g., 
Cr, a simple metal (with a room temperature resistivity [25] 
of about 0.01 mil ■ cm) which orders into a spin-density-wave 
ground state. 
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of the electronic states near the Fermi energy, and define the 
incoherent part through d^ aa = d^ a + d!^° . Therefore, 
unlike the full d-electron operator, d^ a does not satisfy the 
fermion anticommutation rule. Indeed, its spectral function 
integrated over frequency defines w. We therefore introduce 
Ckao- = s0 that Ckan has a total spectral weight 

of 1 and satisfies {cw a a,c^ aa } = 1- 

We then have the effective low-energy Hamiltonian terms 
for the coherent itinerant carriers (H c ) and for their mixing 
with the local moments (H m ): 

He — / electa e^ aa Ck a( j — W y ^ E]^ a( jC^ aa Clicta 

k,CK,<T k,c*,0" 
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Fig. 1. Single-electron spectral function as the sum of co- 
herent and incoherent parts. The single-electron density of states (DOS) 
is plotted against energy (E); Ep 15 the Fermi energy. Each peak may contain addi- 
tional structure due to the multi-orbital nature of the iron pnictides. The percentage 
of the total spectral weight that belongs to the coherent part is defined as w , which 
goes from 1 when the interaction is absent, to when the interaction reaches and 
goes beyond the Mott-transition threshold. 

Effective Hamiltonian 

To study the magnetism, the incoherent spectrum is naturally 
described in terms of localized magnetic moments, leading to 
a matrix J1-J2 model [9]: 

Hj = Jl^S^a ■ Sj,0 + J^Si.a ■ Sj./3 

(ij) 

+ Jh 22 s,, a • s it g. [1] 

Here, Ji and J2 label the superexchange interactions between 
two nearest-neighbor (n.n.,(ij)) and next-nearest-neighbor 
(n.n.n.,((y))) Fe sites, respectively. Both are matrices in the 
orbital basis, a, (3 with these indices summed when repeated. 
Jh is the Hund's coupling. 

Eq. (1) reflects the projection of the full interacting prob- 
lem to the low-energy subspace when the system is a Mott 
insulator (w — 0) and the single-electron excitations have 
only the incoherent part. When the single-electron excita- 
tions also contain the coherent part (w being non-zero but 
small, see Fig. 1), these coherent electronic excitations must 
be included in the low-energy subspace as well. 

We will use the projection procedure of Ref. 26 to con- 
struct the effective low-energy Hamiltonian. We denote by 
d^aa the d-electron operator projected to the coherent part 



Here t labels the three Pauli matrices. In the projection pro- 
cedure leading to Eq. (2), we keep dkoo- as part of the low- 
energy degrees of freedom; the prefactor w in the first equa- 
tion comes from the rescaling Ck Q(T = (l/v^^kalx an d Ek arT is 
therefore the conduction-electron dispersion at w = 1. At the 
same time, we integrate out the high energy states involved in 
Cf- To the leading order in w, this procedure is carried out 
at the w = point which is taken to have a full gap [26]; as a 
result, the effective coupling GkqQ/37 is of order w°. Beyond 
the leading order in w, the coupling constants will acquire 
further corrections. The computation of these corrections is 
difficult, since, at those orders, the spectrum becomes contin- 
uous from the coherent to incoherent part (cf. Fig. 1); it is 
left for future work. Still, our leading-order analysis captures 
the form of the low-energy effective Hamiltonian, which is 

H e f f — Hj + H c + H m . [3] 
Jl-J2 competition 

The superexchange interactions in the iron pnictides contain 
n.n. and n.n.n. terms because of the specific relative locations 
of the ligand As atoms and Fe atoms [9, 16, 27]. To assess 
the tunability of Ji and J2, we consider an over-simplified 
case, illustrated in Fig. 2. Here, only one Fe 3d orbital is 
considered. We assume that the 3d orbital on each of the 
four corners of a square plaquette has an identical hybridiza- 
tion matrix element, V, with one As 4p orbital located above 
the center of the plaquette. The superexchange interaction 
is found to be hj oc X^ r En s ( r )] 2 ' where r labels a pla- 
quette in the 2D square lattice and the summation is 
over the four Fe sites of a plaquette. For classical spins, this 
is the canonical case of magnetic frustration: all states with 
s(r) = are degenerate. Written in the form of Eq. (1), 
this corresponds to J2 = Ji/2. This discussion is instructive 
for the understanding of the realistic exchange interactions in 
the iron pnictides. Several aspects are neglected in the sim- 
plified analysis given above. First, multiple 3d orbitals are 
important, and the hybridization is orbital-sensitive. Both 
the Ji and J2 interactions are therefore matrices. Second, the 
real band structures must be described by more complex d-p, 
p-p, and d-d tight-binding parameters. Both features spoil the 
elementary J2 = Ji/2 relationship. Still, the simple consid- 
erations given above suggest that the overall strength of J2 
and Ji - i.e. the largest eigenvalues of the J2 and Ji matri- 
ces - are comparable with each other. Detailed analysis of 
the matrix elements indicates that there are more entries in 
the J2 matrix than in the Ji matrix that correspond to the 
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Fig. 2. A square plaquette of Fe ions, with an As ion sitting 
above or below the center of the plaquette. The hybridization 

term is written as hv = 5Zr V[pl-( r ) X)n dv( r ) + h.C.], and the energy levels 
for the Fe 3d orbital and As 4p orbital are and e p , respectively. The resulting 
superexchange interaction is hj = [2V 4 /(e p - e d ) 3 } J] r E n s(r)] 2 . 



dominating antiferromagnetic component, and that the over- 
all magnitude (the largest eigenvalue) of the J2 matrix will 
be somewhat larger than half of that of the Ji matrix. This 
conclusion is supported by the fitting of the ab initio results 
of the ground state energies for various magnetic configura- 
tions in terms of Ji and J2 parameters in the (non-matrix) 
Heisenberg form [16, 27]. This range of J2/J1 leads to a two- 
sublattice collinear antiferromagnetic ground state, consistent 
with the results of the neutron scattering experiment [2]. 

On the one hand, the above argument implies that the 
magnetic frustration effect is strong, and can provide signifi- 
cant quantum fluctuations leading to a reduced ordered mo- 
ment. On the other hand, it suggests that the degree to which 
J2 / Ji can be tuned in practice could be limited. 



Magnetic quantum critical point 

The order parameter for the two-sublattice antiferromagnet 
appropriate for J2/J1 > 1/2 is the staggered magnetization, 
m, at wave vector Q = (tt, 0). The effective theory for the Hj 
term alone corresponds to a (f> 4 theory whose action is of the 
form S ~ r<f> 2 + ucf> 4 . The coupling to the coherent quasipar- 
ticles is given by given by the H m term of Eq. (2); it causes a 
shift of the tuning parameter r and also introduces a damping 
term. These contributions to the r coefficient are given by: 



Ar + iF 



2 2 _/(ek+ q , a ) - /(ek,/3) 

i/kqa/37^'7 . 



[4] 



Here /(e) is the Fermi-Dirac distribution function and a 7 is 
an orbital-dependent coefficient: a 7 s 7 appears in the or- 
der parameter for the (tt, 0) antiferromagnet. Note that both 
<?k,qa/37 and ek./3,£k+q,Q are linear order in w. We can infer 
from Eq. (4) that the damping term is of the order w° at low 
energies: for |o> n | <IC wW (W is the bandwidth), V — j\(jJ n \, 
where 7 is, to leading order in w, the constant value associated 
with the couplings and density of states of the w = 1 case. 
Note that 7 is non-zero because, for the parent compounds, 
Q connects the hole pockets near the V point of the Brillouin 
zone (BZ) and the electron pockets near the M points (in the 
unfolded BZ notation). At the same time, 7 does not diverge 
since the nesting is not perfect. The existence of the linear 
in u) damping term is in contrast to the doped case, where Q 
no longer connects the hole- and electron- Fermi surfaces [11]. 
Importantly, we can also infer from Eq. (4) that the leading 
frequency and temperature independent term Ar = wAq is 



linear in w, with A Q = Ek, a ,/3, 7 G kqQ(37 a 7 [e(£F - -Bk+q) - 
Q(Ef — -Eic)]/(-Ek,/3 — -Ek+q.a) (where O is the Heaviside func- 
tion) is independent of w, and positive. 

The low energy Ginzburg-Landau theory then takes the 
form, 

S = J dqjdcu [r(w) + c(q — Q) 2 + lu 2 +j\u>\] [m(q, uj)] 2 

+ u W^JdmJduJisi^^c^ 8^2^ [m] 4 + . . . ,[5] 

where r(w) = r(w — 0) + wAq. r(w — 0) is negative, placing 
the system at w — to be antiferromagnetically ordered. The 
linear in w shift, wAq, causes r(w) to vanish at a w — w c , 
leading to a quantum critical point. In terms of the external 
control parameter S, shown in Fig. 3, w = w c defines S — 8 C . 
The (f> 4 theory describes a z — 2 (where z is the dynamical 
exponent) antiferromagnetic quantum phase transition, which 
is generically second order. 

The 0(3) vector m, corresponding to the (tt, 0) order, is 
accompanied by another 0(3) vector, m' that describes the 
(0, tt) order. These two vector order parameters accommo- 
date a composite scalar, m • m', the order parameter for an 
Ising transition [28, 10, 11]. In turn, the Ginzburg-Landau 
action, Eq. (5), contains a quartic coupling -£t(m • m') 2 [as 
well as «'m 2 (m') 2 ]. In the 2 = 2 case here, the 4> 4 theory is 
at effective dimension d + z — 4. At the QCP of the 0(3) 
transition, the u quartic coupling term is marginally relevant 
in the renormalization group sense. The T — transition 
could therefore either be turned to first order, or be split into 
two continuous transitions, one for the Ising transition, whose 
scalar order parameter is m ■ m' , (which corresponds to the 
structural distortion when it is coupled to some structural de- 
grees of freedom) the other for the 0(3) magnetic one. Either 
effect is expected to be weak, due to the marginal nature of 
the coupling. 

The magnetic quantum criticality will strongly contribute 
to the electronic and magnetic properties in the quantum 
critical regime. We note that since d = z — 2, there are 
(marginal) logarithmic corrections to simple Gaussian critical 
behavior [29]. Following discussions in, e.g., Ref. [29], we ex- 
pect that the specific-heat coefficient will be C/T ~ ln(l/T), 
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Fig. 3. Magnetic quantum phase transition in the parent 
compounds of the iron pnictides. The blue solid/black dashed lines 
represent the magnetic/structural transitions, respectively. 5 is a non-thermal control 
parameter: increasing S enhances the spectral weight in the coherent part of the 
single-electron excitations (cf. Fig. 1). The QCP, at 8 = 8 C , separates a two- 
sublattice collinear AF ground state from a paramagnetic one. A specific example for 
8 is the concentration of P doping for As: a parent iron pnictide with As is an anti- 
ferromagnetic metal, whereas its counterpart with P is non-magnetic; the possibility 
that the latter is superconducting is not shown in the phase diagram. 
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the NMR relaxation rate 1/Ti oc const, and (in the presence 
of disorder scattering that smears the Fermi surface) the re- 
sistivity p oc T. 

Tuning parameter and variation of magnetic order 

The parent materials of the different iron arsenides will have 
different internal pressures and "c/a" ratios, and will corre- 
spondingly have different ratios of the electron-electron inter- 
action to the effective bandwidth. According to our theory, 
the resulting variation of the coherent spectral weight w will, 
in turn, tune the control parameter r in Eq. (5), and the 
ordered moment will change accordingly across the different 
compounds. 

Neutron scattering experiments have indeed found that 
the ordered moment does vary across the parent arsenides. 
The moment associated with Fe-ordering at low tempera- 
tures is about 0.2-0.3^ /Fe in NdOFeAs [30, 31], 0.4^ S /Fe 
in LaOFeAs [2], 0.5^/Fe in PrOFeAs [32, 33], and 0.8-1.0 
Ms/Fe in CeOFeAs [15], BaFe 2 As 2 [34], and SrFe 2 As 2 [36]. 



Comparison with DFT studies 

We have considered the mechanism for quantum fluctuations 
having in mind the proximity to the Mott limit, where the 
instantaneous atomic moment is large (a few /j,b /Fe) to begin 
with. Most DFT calculations have shown that the ordered 
moment in the antiferromagnetic ground state is large, of the 
order 2/us/Fe. Moreover, such a large ordered moment was 
found within DFT not only for the parent iron pnictides, but 
also for their doped counterparts. 

Since DFT calculations neglect quantum fluctuations, we 
are tempted to interpret the large DFT-calculated moment 
as essentially the instantaneous atomic moment. Quantum 
fluctuations will then lead to a reduced ordered moment in 
the true ground state. The Ji-J 2 competition together with 
the coupling of the local moments to the coherent itinerant 
electronic excitations arising naturally in the Mott-proximity 
picture we have described, is just such a mechanism for quan- 
tum fluctuations. 



Asx — <5P<5 Series of the Parent Iron Pnictides 

Since the c-lattice constant in LaOFeP is smaller than that 
in LaOFeAs, these considerations suggest that the coherent- 
electron spectral weight of the iron phosphides is larger than 
that of the iron arsenides. A consequence is that, in con- 
trast to the arsenide, the phosphide does not have a magnetic 
transition. [35] We then propose that a parent iron pnictide 
series created by P doping of As presents a means to unmask 
a magnetic quantum critical point. Our purpose is better 
served the weaker the superconductivity is in the P end ma- 
terial. LaOFeAsi-^Pa is promising, since LaOFeP is a weak 
superconductor whose T c is only a few Kelvin or may even 
vanish [37, 38, 39]. CeOFeAsi_aPi may also be of interest 
in this context. While CeOFeAs [15] is antiferromagnetic, 
CeOFeP is a paramagnetic metal [40]. We remark in passing 
that P-doping for As is more advantageous than external pres- 
sure, because the latter is known to cause a volume collapse 
[41]. It would be interesting to search for a substitution for 
As such that w could be reduced, leading towards to the Mott 
insulating state. 

To understand the tuning of the microscopic electronic 
parameters, we have carried out density-functional-theory 
(DFT) calculations on both CeOFeAs and CeOFeP for com- 
parison. We find that the d-p hybridization matrix is larger in 
CeOFeP than in CeOFeAs. This is consistent with the quali- 
tative consideration that, compared to CeOFeAs, CeOFeP has 
a higher internal pressure and, hence, a higher kinetic energy 
and smaller ratio of the interaction to the bandwidth, thus a 
larger coherent weight w. 



Discussion 

We have developed a framework to describe the quantum mag- 
netism of the iron pnictides, appropriate for electron-electron 
interactions that are of an intermediate strength to place the 
materials at the delicate boundary between itinerancy and 
localization. Our description takes into account the inter- 
play between the itinerant and local-moment aspects, which 
are naturally associated with the interaction-induced coher- 
ent and incoherent parts of the electronic excitations. En- 
hancement of the spectral weight associated with the coher- 
ent electronic excitations weakens the magnetic order, and 
induces a magnetic quantum critical point. Our characteri- 
zation of the magnetic excitations is important not only for 
the understanding of the existing and future experiments in 
the normal state, but also for the microscopic understanding 
of high temperature superconductivity in the iron pnictides 
and related metallic systems close to a Mott transition. In 
addition, realization of a magnetic quantum critical point in 
the iron pnictides provides a new setting to explore some of 
the rich complexities [42, 43] of quantum criticality; this is 
much-needed since quantum critical points have so far been 
explicitly observed only in a very small number of materials. 
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